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We show the existence of massive Dirac fermions in electronic band structures of a few Fe atomic 
layers with perpendicular magnetization. Based on a tight binding model fitted to ab-initio band 
structure, we observe four distinct massive Dirac fermions near the Fermi level, which result from 
atomic spin-orbit coupling of Fe and a band inversion between Fe 4:S-3d„,2_y2 hybrid orbital band 
and 3dxy orbital band. These lead to a valence band with finite Chern integer (-1-2) and chiral edge 
modes near the Fermi level. When the chemical potential is set inside the Dirac gap by carrier 
doping, the Hall conductivity exhibits a plateau-like structure with quantized value 2(e^/h), and 
orbital magnetization shows a prominent increase, latter of which is mostly due to chiral orbital 
motion of electrons along the edge modes. We discuss the stability of the Dirac fermions in Fe(OOl) 
monolayer on MgO(OOl) substrate and Fe(OOl) bilayer case. 

PACS numbers: 


I. INTRODUCTION 

Since the discovery of grapheneji*^ a number of novel 
two-dimensional electronic phases have been theoreti¬ 
cally proposed4i^ Such efforts include a proposal of so- 
called ‘Chern insulator’ in graphene with ferromagnetic 
substrates. Thereby, it was theoretically proposed that 
magnetic proximity effect from the substrate in combi¬ 
nation with atomic spin-orbit interaction of carbon atom 
could give rise to a finite mass gap in the Dirac fermion of 
graphene, resulting in the quantized Hall conductancei^ 
On the experiment end, a few atomic layers of 3d-electron 
ferromagnets themselves such as Fe, Co, Ni and their al¬ 
loys can be easily grown on various substratesfi When 
the thickness of magnetic layers reaches atomic scale, 
magnetic anisotropy energy associated with surface mag¬ 
netism dominates over magnetostatic energy, often mak¬ 
ing magnetic easy axis to be perpendicular to the layer; 
perpendicular magnetic anisotropy (PMA). — To ob¬ 
tain comprehensive understandings of the PMA phenom¬ 
ena, it is also important to extract common feature in 
electronic band structures of a few atomic layers of d- 
electron ferromagnets such as Fe, Co and Ni. 

In this paper, we show the existence of massive Dirac 
fermions near the Fermi level in electronic band struc¬ 
tures of bcc Fe(OOl) a few atomic layers with perpendic¬ 
ular magnetization (Fig.[T|). The Dirac fermions are well 
described by a three-orbital tight-binding model com¬ 
posed of 3dxy, 3dx2-y2 and 4s orbitals of Fe atom, where 
an atomic spin-orbit interaction of Fe and a band inver¬ 
sion between Zdxy orbital band and 4s-3da;2_j,2 hybrid- 
orbital band play essential role for the emergence of the 
Dirac fermions. From ah-initio band calculations, the 
mass of the Dirac fermions is estimated to be around 
80/60 meV for a Fe monolayer without/with MgO(OOl) 
substrate. In a free-standing Fe(OOl) monolayer, the 
massive Dirac fermion results in plateau-like feature in 
Hall conductivity and dMorb/d/i near the Fermi level. 


where Morb out-of-plane orbital magnetization and /r the 
chemical potential. In experiments, the chemical poten¬ 
tial can be controlled by an electric gate voltage (out-of- 
plane electric field). It is shown that the Dirac fermions 
in Fe monolayers are robust against broken out-of-plane 
inversion symmetry induced by the electric voltage. A 
comparison with existing ab-initio band calculation o^^i^^ 
suggests that the Dirac fermions of the same origin can 
be also found in electronic band structures of bcc Co(OOl) 
and hep Co(lll) monolayers. 



FIG. 1: Four massive Dirac fermions in the electronic band 
structure of a perpendicularly magnetized free-standing Fe 
monolayer 
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FIG. 2: (a) Electronic band structures for minority spin near 
the Fermi level (without atomic SOI) from LSDA calcula¬ 
tion (black) and 8-bands tight-binding model (blue). The 
Dirac fermion is depicted by a green circle, (b) from the 
6-bands tight-binding model (with atomic SOI). The low¬ 
est band which acquires the Chern number 2 is drawn with 
red bold line. The massive Dirac fermions are depicted by 
a green circle, (c) from LSDA calculation with MgO(OOl) 
substrate and without atomic SOI. The orbital characters of 
bands labeled as 1,2 are 3dxz, Sdyz, while those of 3, 4, 5, 6 
are 3dj,2_j,2, Sd^y, 4s, Sd ^2 respectively. 


II. ELECTRONIC BAND STRUCTURE 

Ab-initio calculations are performed using film full- 
potential linearized augmented plane wave (FLAPW) 
method^S^^ based on the local spin density approxi¬ 
mation (LSDA),— in which the core states are treated 
fully relativistically and the valence states are treated 
semi-relativistically. LAPW functions with a cutoff of 
|k-|- G| < 3.9 a.u.“^ and muffin-tin sphere radii of 2.2, 
2.2 and 1.4 a.u. for Fe, Mg and O atoms are used, 
where the angular momentum expansion inside the MT 
spheres is truncated at ^ = 8 (Fe and Mg) and i = Q 
(O) for the wave functions, charge density and poten¬ 
tial. The Fe/MgO was modeled by an Fe monolayer on 
a six-atomic-layer MgO(OOl) substrate, where Fe atoms 
locate on top of the O atoms, assuming the in-plane lat¬ 
tice constant matching to the calculated value of bulk 
MgO while the out-of-plane coordinates of Fe atoms are 
fully optimized using the atomic-force FLAPW calcu¬ 
lations. Note that the structural parameter obtained 
within LSDA qualitatively agree with those analyzed by 
surface X-ray diffractioit^^ as demonstrated previously;^ 

The ab-initio electronic band structure for minority- 
spin bands of free-standing Fe monolayer with perpen¬ 
dicular magnetization is fitted near the Fermi level with 
a tight binding model composed of five 3d orbitals, 4s, 
4px and Apy orbitals of Fe atom; 

i,m {{i,j))\mn 

+ Ago 

i\mn 

( 1 ) 

mt (D,mt) is the electron creation (annihilation) oper¬ 
ator on site I, for orbital m, and spin up (where spin- 
quantization axis is taken along the perpendicular mag¬ 
netization direction, z-direction). The first term repre¬ 
sents an effective atomic energy of each orbital, which 
includes all the on-site type energies felt by respective 
orbital, such as exchange splitting energy from major¬ 
ity spin electrons and crystal fields caused by surround¬ 
ing electrons. The second term represents neighboring 
hopping term with inter-atomic/intra-atomic transfer in¬ 
tegrals tij-rnnt- The angular bracket in ((*,j)) stands 
for summation of nearest neighboring and next nearest 
neighboring sites. The third term is the atomic spin-orbit 
coupling with coupling strength taken to be Ago = 50 
meV from several literatures. Since an exchange splitting 
between majority spin bands and minority spin bands es¬ 
timated from ab-initio band calculations is about 3 eV, 
being 60 times larger than the atomic spin-orbit cou¬ 
pling strength of Fe, we consider only diagonal part of 
spin-orbit interactions with respect to spin index; an ex¬ 
pected correction due to the off-diagonal parts, [L±][<S';p], 
is evaluated to be on the order of 1 meV, being fairly 
negligible compared to that from the diagonal part. The 
interatomic transfer integrals tij-mn-f are given by matrix 
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elements in the Slater Koster table;^ such as Vsp, 
Vddai VddTT, VddS and SO On. These matrix elements in the 
Slater Koster table along with the effective atomic ener¬ 
gies Cmt are used as fitting parameters, whose best fitted 
values in the absence of spin-orbit interaction are shown 
in Table. I. The table shows reasonable fitting values 
compared to the Solid State Tabled and the lattice con¬ 
stant of the square-lattice Fe monolayer evaluated from 
the same first principle calculation (a = 2.95A). 

Figure [2ja) plots the band structure along high sym¬ 
metry k points obtained from both our eight-bands tight 
binding model and LSDA calculations with Ago = 0. No¬ 
tably, there exists four linearly dispersive band crossings 
along the zone boundary connecting M point and X 
point near the Fermi level. As shown below, these linear 
dispersions are well described by massless Dirac fermions 
without the spin-orbit interaction. With the atomic spin- 
orbit interaction, each massless Dirac fermion acquires 
same sign of the mass, endowing a valence band with 
finite Chern number +2. 

The four Dirac cones are enveloped by other two dis¬ 
persive bands composed of dxz and dyz, while the out- 
of-plane mirror symmetry allows us to treat these two 
bands separately from the other six bands composed of 
s,dx^-y^, dxyi Px: Py and dz2, because {z odd\H\z even) = 
{z odd\azdzHazd'z\z even) = —(z odd\H\z even) = 0. 
The presence of substrate, such as MgO(OOl), breaks 
the out-of-plane mirror symmetry, giving rise to a fi¬ 
nite mixing between these two groups of bands, while, 
more importantly, the oxide substrate endows with rel¬ 
atively strong charging energies the out-of-plane d or¬ 
bitals, such as dxz, dyz and (i 22 .Thus, dxz-dyz bands are 
brought into a higher energy region, while the four Dirac 
fermions composed by 4s, 3dx^-y^ and Sdxy remain in¬ 
tact (FigISlJc)). For the sake of clarity, we consider in the 
following a free standing Fe monolayer and treat these 
two groups of bands separately, unless dictated otherwise. 
The case with MgO(OOl) substrate will be discussed later 
more carefully. 


III. BAND INVERSION MECHANISM 

The emergence of the Dirac fermions and the valence 
band acquiring the Chern number +2 result from (i) a 
band inversion between ^dxy orbital band and dsSdx^-y^ 
hybrid-orbital band, and (ii) a complex-valued band mix¬ 
ing between these two bands mediated by the atomic spin 
orbital interaction. Due to the orbital symmetry, the 
nearest neighbor intra-orbital transfer integrals make the 
hybrid orbital band to have positive curvature at the F 
point and negative curvature at the M point, while the 
2idxy orbital band to have negative curvature at the F 
point and positive curvature at the M point. As a re¬ 
sult, the hybrid-orbital band comes lower than 'idxy or¬ 
bital band at the F point, while comes higher than Zdxy 
at the M point (’band inversion’). 

In the presence of the atomic spin-orbit interaction 


TABLE I: (upper) Tight-binding hopping parameters for mi¬ 
nority spin bands, obtained from the fitting to the LSDA 
calculation. The energy unit is eV. (lower) On-site effective 
atomic energies for minority-spin 3d electrons and 4s, 

4pj| electrons. 



(J 

TT 

<5 

Fe d-Fe d (NN) 

-0.363 

0.261 

-0.061 

Fe d-Fe d (NNN) 

-0.080 

0.064 

-0.025 

Fe p-Fe d (NN) 

-0.392 

0.157 

— 

Fe p-Fe d (NNN) 

-0.118 

0.047 

— 

Fe p-Fe p (NN) 

2.905 

-1.063 

— 

Fe s-Fe d (NN) 

-0.392 

— 

— 

Fe s-Fe d (NNN) 

-0.118 

— 

— 

Fe s-Fe s (NN) 

-1.419 

— 

— 

Fe s-Fe p (NN) 

2.208 

— 

— 


s p 

dxz 

dyz d^2 y2 

dxy 

ei,t 2.315 8.915 

0.350 0.233 

0.234 0.160 

0.463 


with the out-of-plane ferromagnetic moment, dxy orbital 
is mixed with dx 2 -y 2 orbital character with pure imagi¬ 
nary coefficient, i.e. dxy —f dxy + iadx^-y^, with a pro¬ 
portional to the spin-orbit interaction. Because of this 
mixing, the hybrid orbital band and 3dxy orbital band 
acquire a complex-valued band mixing, which takes form 
of t' sin kx sin ky + iat(cos kx — cos ky) in the momentum 
space; sin kx sin ky and coskx — cosky are from the 3dxy 
and 3dx2-y2 orbital symmetry respectively. As shown 
below, the band inversion and the complex-valued band 
mixing result in the valence band with Chern number 
-\- 2 . 

The valence band with Chern number +2 can be well 
described by the lowest energy band of a three-bands 
tight-binding model composed of 4s, 3dxy and 3dx2-y2 
orbitals. In the momentum space, the model takes form 
of. 




( Es{k) -2y/W!,y^SxSy VSVsdicx - Cy) 

* Exy(^k) 2i\so 

* * Ex2_y2{k) 


( 2 ) 


where E^^k) Si -t- tii^Cx Cy) F Cx^y — coskx^y, 

Sx,y = sinkx,y, subindex i = s,xy,x^ — representing 
4s, 3dxy, 3dx2-y2 orbitals. Si being their effective atomic 
energies, ti and t) denotes their nearest and next nearest 
neighbor intra-orbital transfer integral. The lowest band 
energy of eq. ([H will be denoted as (fc) henceforth. 

To see the band topology for the lowest band, let us 
first derive an effective 2x2 Hamiltonian out of H^^^{k). 
To this end, notice first that the lowest band energy is 
always smaller than Es{k), Exy{k) and Ej. 2 _y 2 {k) for all 
fc; E_ < nnn{Eg, Exy, Ex 2 _y 2 ). In such an occasion, we 
may describe the lowest band, by treating either one of 
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these three orbitals as a high energy degree and deriving 
an effective 2x2 Hamiltonian for the other two; the low¬ 
est band of is identified with that of the 2 x 2 

Hamiltonian. As shown in Fig. [2Ib) and Fig. [3J the fo¬ 
cused valence band has mainly 4s and 3dxy characters. 
Thus, we regard orbital as the high energy degree 

of freedom and treat its couplings with 4s and 3dxy as 
perturbations. This leads to the 2x2 effective Hamilto¬ 
nian for ‘4s band’ and 5dxy band. Quantitatively speak¬ 
ing, ‘4s band’ thus introduced is a 4 s-3dx2-y2 hybrid- 
orbital band rather than purely 4s orbital band, because 
of larger mixing between 4s and 3 dx2-y2 coming from 
Vsd Ago)- In fact, the hybrid band has larger 8 ^ 3 , 2 _j ,2 
orbital character than 4s orbital in a certain momentum 
region (see Fig. [3|). 

The degenerate perturbation theory gives the 2x2 
Hamiltonian as follows. 




E Ex2 — y2(l?dj 


with i, j = s, xy and 


Hn = 


( Es{k) -2V3V'^SxSy 

-2V3V'^SxSy Exy{k) 

\ 0 0 



/ 


H, = 


y '\/3\^sd(^^x 2zA 

Equivalently, 


0 V^^sdi^^X Cy) 

0 -2zAso 
0 


(4) 

(5) 


with 


tt2x2 _ I Es{k) Es^xy{k) \ 

“ I X^xy{k) Exy[k) i ’ 


■p — P _ Cy) — _ p _ 4 Aso 

I^S — -^s J^xy—^xy ’ 


( 6 ) 


Eb.xv = —2-\/3y( 


i2V3X SO ^sd{Cx Cy^ 


sd^oc^y 


AE 


The perturbation treatment is valid as far as AE = 
Ex2-y2{k) — E \B positive. This condition is satisfied 
for any k when E = E_{k). 

The 2x2 Hamiltonian has two eigenvalues whose 
smaller one corresponds to the lowest band energy of 
^ 3 x 3 (^)_ Thereby, E-{k) can be obtained from the fol¬ 
lowing self-consistent equation of E -; 


2E_{k) — Es{k) -\- Exy{k) 

- ^(EB{k) -Exyik))^ + 4\EB,xyikW. 


where the right hand side is given by E_{k) itself by way 
of AE = Ex 2 _y 2 {k) — E_{k). The equation has a solution 
for E_ which always satisfies E_ < min(Es, Exy, Ex2-y2) 



FIG. 3: ‘Fat-band’ picture obtained from the 6-bands tight- 
binding model. Respective orbital character is depicted by 
the line width (red) in (a) for 4s, in (b) for 3dj.2_y2 and in 
(c) for 3dxy. 


for any fc, justifying a posteriori the validity of the per¬ 
turbative treatment. 

With this justification in mind, we can readily identify 
the band topology of the lowest energy band of H^^'^{k) 
as that of H^^^{k). As is clear from Fig. [2Kb), the hybrid 
band comes lower than 3dxy band at the F point; Es < 
Exy at fc = F, while otherwise at the M point; Eg > Exy 
dX k = M . Between F and M point, these two bands 
have a mixing due to a finite off-diagonal matrix element 
Es,xy{k). Importantly, the matrix element has a complex 
phase, which acquires 47 r phase, whenever fc goes around 
the F point (or M point). 

f Vk{a.rgEs,xy{k)} ■ dk = An, (7) 

Jas 

where dS denotes an arbitrary loop which encompasses 
the F-point (or M-point). 

The 47 r phase winding of the inter-band matrix ele¬ 
ment and the band inversion between the hybrid band 
and 3dxy band endow the lowest band with the Chern 
number -|-2. To see this, expand the 2x2 Hamiltonian in 
terms of the Pauli matrix, E[^g'^{k) = ao{k)ao + o,{k)-cr^ 
from which the normalized vector n(fc) is introduced by 
n = a/|a|. According to the projective representation 
of the Chern invariants'^— the Chern number for the 
lowest band (Ch_) is given by an integral of solid angle 
subtended by the unit vector over the first Brillouin zone; 
Ch_ = ^ dkxdkyn{k)-{dk^n{k)xdkyn{k)). The in¬ 
tegral is quantized to be integer, which counts how many 
times the unit vector wraps the unit sphere when the mo¬ 
mentum fc wraps the first Brillouin zone once. Now that 
Eg < Exy and Eg^xy = 0 at the F point while Eg > Exy 
and Eg^xy = 0 at the M point, the unit vector points to 
the south pole/north pole of the unit sphere when fc at 
the T/M point respectively. On the one hand, eq. ([7]) 
means that the unit vector always winds twice around 
the pole when fc rotates once around the F point. This 
dictates that the Chern integer for the lowest band is +2. 

The lowest band forms four distinct Dirac fermions 
along the Brillouin zone boundary, {iTi, K 2 , K 4 } = 

{(tt, A"), (tt, — AT), (AT, tt), (—AT, tt)}. Around each Dirac 
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FIG. 4: Schematic picture of dual magnetic magnetic charges 
and dual magnetic field in the three-dimensional space sub¬ 
tended by kx, ky and Aso- The Chern integer for each Aso is 
given by the surface integral of the field over the first Brillouin 
zone for the constant Aso-; 


point, the effective 3x3 Hamiltonian conceiving Dirac 
fermion is linearly expanded in small Qy and Ago, e.g. 

H^g^{k) = Mi{qxa)(Tx + M2Xso(Ty + M3{qya)cr^ H-(8) 

with k = q + Ki, and Mi = —0.34eV, M2 = —0.97, 
Mg = —0.66eV for the tight-binding parameters in Table 

I. 

The four Dirac points play role of dual magnetic 
monopoles in a three-dimensional parameter space sub¬ 
tended by kx,ky and Ago'^^— The corresponding 
magnetic field B_(fc,Ago) is associated with a Bloch 
wavefunction for the lowest band, |M_(fc,Ago)) with 
H{k, Xso)\u-) = E-\u-). The magnetic field is a 
rotation of a three-component gauge field A_(fc,Ago), 
J5_(fc,Aso) = V X A_(fc,Aso) with V = {dkx,dky,dxj). 
A_{k, Ago) are gauge connections of the Bloch wavefunc¬ 
tion: A_ = i(M_|V|M_). Due to the four-fold rotational 
symmetry, dual magnetic charges at four Dirac points 
have the same quantized strength 2tt, where their sign is 
same as —sgn[MiM2M3] (Fig. 0]), 

4 

V • B^{k, Ago) = 27r^ S{Xso) 6 {k - K,), (9) 

i=i 

The Chern integer for the lowest band is the total mag¬ 
netic flux penetrating through the constant Ago plane in 
the 3D space, 

Ch_(Aso) = J ^(dkxA-^y{k,Xso) - dkyA_^a;{k,Xso)y 

When Ago goes across the Ago = 0 plane, the Chern inte¬ 
ger changes by —4 (Fig. |4]), 



FIG. 5: (upper/lower) Hall conductivity/orbital magneti¬ 
zation as a function of the chemical potential for 8-bands 
tight-binding model for a free-standing Fe monolayer. Ow¬ 
ing to the out-of-plane mirror symmetry, the contribution can 
be decomposed into that from the 6-bands electronic states 
(Sdxy, Sdx. 2 _y 2 , 3d^2, 4s, red color) and that from 

the 2-bands electronic states {3dxz, Sdyz'. blue color). The 
direct band gap region associated with the massive Dirac 
fermions are specified by a grey-hatched energy window, 
[-0.05 eV,0.03eV]. 


IV. HALL CONDUCTIVITY AND ORBITAL 
MAGNETIZATION 


The hallmark of the existence of massive Dirac 
fermions is the transverse conductivity i^^i^^'^^ The Hall 
conductivity as a function of the chemical potential is 
calculated for the eight-bands tight binding model for a 
free standing Fe monolayer (Fig. [5|) with; 


Ch_(Ag„<0)-Ch_(Ago>0) = -4. (10) 

The time reversal symmetry connects the spinless tight- 
binding Hamiltonian for Ago > 0 and that for Ago < 0 
with the relation i7*(fc,Ago) = i7(— fe, — Ago), which 
leads to Ch_(Ago) = —Ch_(— Ago). Combing this with 
eq. m, we have Ch_(Aso > 0) = +2. 



where = i{un\dky\un) with i7(fc)|u„(fc)) = 

£n\un{k))^ n the band index. Due to the out-of-plane 
mirror symmetry, the conductivity can be decomposed 
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FIG. 6: (a) Electronic band structure obtained from the tight- 
binding model with a finite slab geometry and periodic bound¬ 
ary condition along a;-direction. The black colored points are 
extended over the two-dimensional bulk, while the red/blue 
color points are for those eigenstates localized at one/the 
other boundary {y = 0 ot y = L) with L = 200. Two chiral 
edge modes are terminated by a bulk band near k = 0. (b) 
Schematic picture of an expected band structure (b) without 
the bulk band near k — 0 (c) with the bulk band near k = 0, 
where two chiral edge modes (partially) go across the direct 
band gap associated with massive Dirac fermions. 


into the 2-bands contribution (from 3dzx and 3dyz or¬ 
bitals) and the 6-bands contribution (from 4s, 3dx2-y2, 
3dxy, 3(iz2, Apx, Apy orbitals). When the chemical po¬ 
tential is inside the Dirac gap (^ ~ 0), the Hall conduc¬ 
tivity shows a prominent peak structure with a maxi¬ 
mum value around 2e^//i. The peak structure is mainly 
due to a nearly quantized contribution from the 6-bands 
electronic states. The quantized value is approximately 
2 e^//i which is a direct consequence of the four mas¬ 
sive Dirac fermions near p = 0. A slight deviation from 
the quantization is attributed to another small but non¬ 
vanishing dual magnetic fields associated with a disper¬ 
sive band near the F point. 

The 2-bands electronic state also gives a nearly quan¬ 
tized contribution jh to the Hall conductivity near the 
Fermi level (/i ~ —O.leV); 3dxz and 3dyz orbital bands 
comprise another SOTinduced direct band gap at the M 
point. The tight-binding Hamiltonian for 3dxz and 3dyz 
orbitals is expanded linearly in small qx%, 9^ ~ 9^ and 


Ago with fc = q -t- (tt, tt); 

H(k) = Xso (Ty + (YddS ~ ^ddTr){dxO-){qya) CTx 


where and denote the Slater-Koster hopping 
parameters between the (next) nearest neighboring Fe d 
orbitals. The expansion dictates that the dual magnetic 
fields for the two bands have dvr magnetic charge at the M 
point on the Xso = 0 plane. This in combination with the 
symmetry property H*{k,Xso) = H{—k,—Xso) requires 
that the integral of the dual magnetic field {dk^An^y — 
dkyAn,x', n is either lower or higher band out of the two 
bands) near the M-point is quantized to be 27r in the 
smaller Xso limit. When the chemical potential is inside 
the SOI-induced gap at the M point, one of the two bands 
is partially filled while the other is empty. Since a fermi 
surface associated with the hlled band is large enough 
compared to a distribution of the magnetic flux around 
the M point, the Hall conductivity from the two-band 
electronic state is nearly quantized to be e^/h as in Fig.[S] 

Emergence of the massive Dirac fermions also results 
in peculiar chiral modes localized near the boundary of 
a two-dimensional Fe monolayer. Fig. [6] shows an elec¬ 
tronic band structure of the 6-bands tight-binding model 
with periodic/open boundary condition along the xjy- 
direction of the square-lattice Fe monolayer. When pro¬ 
jected onto a surface crystal momentum axis, the four 
massive Dirac fermions at fc = Ki, K 2 , K 4 ^ reduce 

to three distinct valleys with a direct band gap, located 
at kx = TT, ±K respectively. Now that the gap endows the 
lower bulk band with the Chern number +2 as described 
above, the bulk-edge correspondence^^— dictates that 
two localized chiral edge modes appear in the direct band 
gap of the three valleys (Fig.lHb)- In the present case, the 
direct band gap is also masked by another dispersive bulk 
band located at F-point, mainly composed of 3dx2-y2 or¬ 
bital (Fig. [2|). As a result, the chiral edge modes are 
terminated by the dispersive bulk band around kx = 0 
(Fig. [nii,c). 

The chiral modes give rise to large out-of-plane orbital 
magnetization when the chemical potential p set inside 
the Dirac gap. When increasing p inside the gap, elec¬ 
trons are added up into the edge modes, which enhances 
chiral electric currents flowing around the boundary of 
the two-dimensional system. Irrespective of details of 
the energy dispersion of the chiral modes, the increase 
of the current is proportional to the increase of p. Such 
chiral edge current contributes to a macroscopic orbital 
moment (r x p)^r^ which results in a linear increase 
of the magnetic moment with respect to the carrier dop¬ 
ing near the Fermi level. To see this situation, we have 
calculated the orbital magnetization based on the Streda 
formula—*^ The magnetization is the derivative of the 
free energy in the magnetic field H, while the total num¬ 
ber of electrons N is the derivative in p. This leads to 
dN/dH = dM/dp, provided that the free energy is ana¬ 
lytic in p and H. According to Streda;^ dN/dH can be 
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FIG. 7: (a) Electronic band structures for minority spin near 
the Fermi level for a free-standing Fe bilayer with atomic 
SOI (black) and without atomic SOI (red), obtained from 
tight-binding calculations. The Dirac fermions are depicted 
by green dotted circles. (b)(c) Hall conductivity and orbital 
magnetization as a function the chemical potential /r. (b) fj, 
ranges from -0.4 eV to -0.2 eV (c) /r ranges from 0.1 eV to 0.3 
eV. 


expressed only in terms of the current operators; 


dM 

dfj, 


1 

ec 


- -TT[J,G+iii)JySifi -%)- h.c.] I (11) 


where Ji, being the current operator {v = x, y) and 
G^(/r) = l/[(/r ± i5)I — V] and % being lesser, greater 
single-particle green functions and Hamiltonian respec¬ 
tively. By an integration over y, the orbital magneti¬ 
zation is calculated from the eight-bands tight binding 
model for a free standing Fe monolayer (Fig. [5]). Like 
the Hall conductivity, the result is decomposed into the 
2-bands and the 6-bands contributions. The calculated 
magnetization exhibits a significant increase as a function 
of /i when y, is set inside the Dirac mass gap. The break¬ 
down into the two contributions shows that the increase 
is mainly due to the 6-bands electronic states, indicat¬ 
ing that the orbital moment near y = 0 mainly comes 
from an orbital motion of electrons along the chiral edge 
modes. In fact, dM/dy is nearly quantized in the unit 
of e/hc, which counts the number of chiral edge modes 
inside the Dirac gap. 


V. SUBSTRATE, ELECTRIC GATE VOLTAGE 
AND MULTIPLE-LAYER EFFECTS, AND HCP 
CO(lll) MONOLAYER CASE 

The Dirac fermions found in a free-standing Fe mono- 
layer are robust against various perturbations such as 


the oxide substrate, out-of-plane applied electric field 
(e.g. electric gate voltage applied perpendicular to the 
layer) and multiple-layer effects. Firstly, being a dou¬ 
bly degenerate point in an electronic energy band struc¬ 
ture, the dual magnetic monopole (charge) discussed 
above is a stable point defect in the 3-dimensional pa¬ 
rameter space subtended by k^, ky and Ago; they can¬ 
not disappear by themselves. The symmetry property 
H*{k,Xso) = H{—k,—Xso) further requires these defects 
to be in the Ago = 0 plane, which guarantees the existence 
of massive Dirac fermions even for small Ago. To annihi¬ 
late these Dirac fermions, one generally needs to either 
re-invert the band inversion between 3dxy band and 4s- 
3dx2_y2 hybrid band or reduce completely the inter-layer 
couplings among 4s orbital, idxy and idx 2 _y 2 orbitals. 
Unlike out-of-plane 3d orbitals, however, these in-plane 
3d orbitals has little influences from the substrate and 
out-of-plane electric field. As a result, we can readily 
find the massive Dirac fermions of the same origin even 
in the presence of various perturbations. 

Fig. [5] shows an electronic band structure for Fe (001) 
monolayer with MgO(OOl) substrate, where every Fe 
atom locates right above the oxygen of the MgO sub¬ 
strate. Due to crystal fields from these oxygens, three 
out-of-plane 3d orbital bands are brought up into a higher 
energy region. Due to the charge neutrality, the Dirac 
fermions formed by dxy orbital and 4s-3da;2_j,2 hybrid 
orbitals come lower than the Fermi level {E ~ —0.6eV). 
The size of the SOC-induced Dirac gap is estimated 
around 60 meV from the ab-initio band calculation. The 
applied out-of-plane electric field has little effects on 
these Dirac fermions either. Even under a very large 
out-of-plane electric field (±1 V/A), four Dirac fermions 
are barely affected 

Fig. [7] shows an electronic band structure obtained 
from a tight-binding model for a free-standing Fe bilayer, 
where the number of massive Dirac fermions are doubled. 
Due to interlayer hoppings, a Dirac fermion from one 
layer and that from the other repel with each other in en¬ 
ergy. When the chemical potential is around these Dirac 
gaps, the transverse conductivity shows a peak structure 
with its maximum value around 2e^//i. The out-of-plane 
orbital magnetization increases as a function of the chem¬ 
ical potential inside the gaps. These features are essen¬ 
tially same as in the free-standing Fe monolayer case. 

A comparison between an existing ah-initio band cal- 
culatioiii^ and tight-binding analysis indicates that the 
massive Dirac fermions of the same kinds are also in¬ 
duced by the atomic spin-orbit interaction in minority- 
spin band in a hep Co(lll) monolayer with perpendic¬ 
ular magnetization. Thereby, dxz and dyz orbital bands 
comprise two massive Dirac fermions with positive mass 
at K and K' point respectively, which correspond to 47r 
magnetic charge at the M-point in the Fe(OOl) mono- 
layer case. Meanwhile 4s, dx^-y^ and dxy orbitals form 
six massive Dirac fermions with positive mass along the 
high symmetric k lines connecting F and K{K') and 2 
massive Dirac fermions with negative mass at K and K' 











respectively. These results in a valence band with the 
Chern number +2, where the band has 4s character at 
r and dxy^x‘^-y ‘2 characters at the Brillouin zone bound¬ 
ary. Due to the difference between nominal valence of 
iron atom and cobalt atom, the Dirac fermions in the 
Co(lll) monolayer appear in a lower energy region than 
those in the Fe(OOl) monolayer case. 

VI. CONCLUSION 

Massive Dirac fermions are discovered near the Fermi 
level of an electronic band structure of Fe ultra-thin 
film. The Dirac gap is induced by atomic spin-orbit 
coupling on the order of 50 meV. The topological gap 
opening results from a band inversion between idixy or¬ 
bital band and 4s-3da;2_y2 hybrid orbital band, giving 
rise to a finite Chern number in a valence band. Inside 
the gap, the Hall conductivity (v.s. chemical potential) 


exhibits plateau-like structure with nearly quantized val¬ 
ues, while orbital magnetization (v.s. chemical potential) 
increases rapidly due to the macroscopic orbital moment 
induced by topological chiral edge modes. The magni¬ 
tude of the calculated orbital magnetization due to the 
chiral edge current is on the same order of experimen¬ 
tal literature valued being hardly negligible in general. 
Massive Dirac fermions in Fe ultra-thin film are shown 
to be robust against perpendicular inversion symmetry 
breaking (such as substrates or electric gate voltages) as 
well as multi-layer effect. More importantly, we found 
in Fe(001)/Mg0(001) that the massive Dirac fermions 
are nicely separated from other dispersive bulk bands in 
energy, which may give a useful hint to explore possi¬ 
ble ‘Chern insulator’ in transition metal ferromagnetic 
thin film. Considering richness of the transition metals 
with various substrates, we anticipate that Dirac fermion 
physics and Chern insulators may be observed experi¬ 
mentally in transition metal thin films in future. 
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